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Abstract: The smallness of the neutrino masses may be related to inflation. The 
minimal supersymmetric Standard Model (MSSM) with small neutrino Yukawa cou- 
plings already has all the necessary ingredients for a successful inflation. In this 
model the inflaton is a gauge-invariant combination of the right-handed sneutrino, 
the slepton, and the Higgs field, which generate a flat direction suitable for inflation if 
the Yukawa coupling is small enough. In a class of models, the observed microwave 
background anisotropy and the tilted power spectrum are related to the neutrino 
masses. 



One of the least satisfying features of inflation is the need to introduce the 
inflaton as an ad hoc ingredient. However, some well- motivated theories beyond 
the Standard Model may already have the gauge-invariant scalar degrees of freedom 
necessary for inflation. One recent example is the use of a flat direction in the minimal 
supersymmetric Standard Model (MSSM) as an inflaton jj]]. The model illustrates 
a successful inflation at low scales whose ingredients are testable at Large hadronic 
Collider (LHC) 0] and CMB ||. There are a few other examples of using some 
gauge-invariant degrees of freedom already present in the theory for the purposes of 
inflation g|. 

In any model, the inflaton potential must be very flat, which is suggestive of 
either a symmetry or a small coupling, or both. We will see that this property of 
the inflaton may be related to the smallness of neutrino masses. Identifying such a 
connection could have important ramifications and could lead to a more fundamental 
theory. The model we will use as an example will contain nothing but the MSSM 
and the right-handed neutrinos. We will show that a viable inflation in this model 
favors the correct scale for the neutrino masses. 

Let us consider the minimal supersymmetric standard model (MSSM) with three 
additional fields, namely the right-handed (RH) neutrino supermultiplets. The rele- 
vant part of the superpotential is 

W = W MS su + hNH u L. (1) 

Here N, L and H u are superfields containing the RH neutrinos, left-handed (LH) 
leptons and the Higgs which gives mass to the up-type quarks, respectively. For con- 
ciseness we have omitted the generation indices. We note that the RH (s) neutrinos 
are singlets under the standard model (SM) gauge group. However in many exten- 
sions of the SM they can transform non-trivially under the action of a larger gauge 
group. The simplest example is extending the SM gauge group to SU (3)c x SU (2)w x 
U(1)y x [7(1) b _l, which is a subgroup of 5*0(10). Here B and L denote the baryon 
and lepton numbers, respectively. This is the model we consider here. In particular, 
the Z7(1)b_l prohibits the RH Majorana masses. 

There is a broad range of phenomenologically acceptable scales Ab-l & t which 
the £/(1)b-l can be broken. In the low-energy effective theory this breaking can 
introduce the RH Majorana mass for neutrinos W D yNN, where M oc Ab-l- 

The active neutrino masses that arise from this are given by the usual seesaw 
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relation h 2 (H u ) 2 /M §, §, where (H u ) is the Higgs vacuum expectation value (VEV). 
Although the seesaw mechanism allows for small neutrino masses in the presence of 
large Yukawa couplings, it does not require the Yukawa couplings to be of order 
one. It still allows one to choose between the large Yukawa couplings and large 
Majorana masses on the one hand, and the small Yukawas, small Majorana masses, 
on the other hand. Viable models for neutrino mass matrices have been constructed 
in both limits, including the low-scale seesaw models 0, 13], in which the Yukawa 
couplings are typically of the order of 

h ~ 1(T 12 , (2) 

or the same order of magnitude, as it would have in the case of Dirac neutrinos in or- 
der to explain the mass scale ~ 0(0.1 eV) corresponding to the atmospheric neutrino 
oscillations detected by Super-Kamiokande experiment. As shown in Ref. such 
models can explain all the observed data on the neutrino masses and mixing, and, in 
addition, one can have dark matter in the form of sterile neutrinos with mass of the 
order of several keV ||. Moreover, in some range of parameters, one can also explain 
the observed velocities of pulsars |jTD[ and the baryon asymmetry of the universe [[1^ . 
In addition, the x-ray background produced by decays of the relic sterile neutrinos 



can play a role in the formation of the first stars that lead to the reionization [11 



In view of the viability of the small- Yukawa scenario for the neutrino masses, 
let us consider the flat directions associated with such small Yukawa couplings in 
connection with cosmological inflation. 

The scalar potential of the MSSM has numerous flat directions made up of 
squarks, sleptons, and the Higgs fields [[13] (for a review see |14|). These directions 
are classified by monomials which are invariant under the SM gauge group. Those 
monomials with B — L = will be also D-flat under £/(1)b-Lj while those with 
B-I / must be multiplied by an appropriate number of N superfields. In 
particular, NH U L is now a D-flat direction. We note that, in the absence of a gauged 
£/(1)b-l, there would be two independent flat directions: N and H U L. These flat 
directions could develop different VEVs, could start oscillating at different times, 
and would decay independently from each other. However, when the SM gauge 
symmetry is extended to include a £/(1)b-l, the true D-flat direction is NH n L. 
As we shall see, this is crucial for the proposed model of inflation and a successful 
reheating. Such a gauged £/(1)b-l which is broken just above TeV is compatible 
with all phenomenological constraints. 
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Let us now work in the basis where neutrino masses are diagonalized. There is 
a flat direction NsH u L3 spanned by the VEV of the lower and upper weak isospin 
components of H u and L3, respectively. The scalar field corresponding to the flat 
direction is denoted by 



where the superscripts refer to the weak isospin components. One must now include 
the soft SUSY breaking terms, such as the mass terms and the A-term. The A-terms 
are known to play an important role in Affleck-Dine baryogenesis ||13| , as well as in 
the inflation models based on supersymmetry [I, [RJ . 
The potential along the flat direction is found to be 

TO 2 /i 2 Ah 

V{<P) = + ^ + ^cos (9 + 9 h + 6 A ) 3 , (4) 

where the flat direction mass is given in terms of the soft masses of N3, H u , and L 3 : 
= (jrift; + m H u + m I 3 ) Here we have used the radial and angular components 
of the flat direction 0#+i0/ = \/2</> exp (i9) , and Oh, 6a are the phases of the Yukawa 
coupling h 3 and the A-term, respectively. We note that the above potential does not 
contain any non-renormalizable term at all. 

The last term on the right-hand side of eq. is minimized when cos (8 + 9h + &a) 
— 1. Along this direction, V(<j>) has the global minimum at = and a local mini- 
mum at (p ~ m^/hs, as long as 

4m^ < A < 3V2 m(t> . (5) 

If A > '&\p2m§, the minimum at </>o will become global. However, this will not be 
relevant for our discussion below. As we will point out later, successful inflation is 
only possible if A is very close to 4m«^, which ensures 0o being a local minimum. At 
the local minimum the curvature of the potential along the radial direction is ~ 
and the curvature along the angular direction is positive. Near <ft = the potential 
reduces to: V ~ m^/hf. If is trapped in the false vacuum and its potential energy 
V dominates the total energy density of the Universe, then inflation ensues. The 
Hubble expansion rate during inflation is given by 

H mf (6) 

We note that H in f <C as <fio <C Mp. This implies that, around the local minimum, 
the potential barrier is too high for the flat direction, 0, to jump over. The situation 
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Figure 1: The inflaton potential. The potential is flat near the saddle point where inflation 
occurs. 



is essentially the same as that in the old inflation scenario jT7|, which could yield 
expansion by many e-foldings, but which suffered from the graceful exit problem. 

However, the barrier disappears when the inequality in eq. (]5|) is saturated, i.e., 
when 

A = 4m0 . (7) 

Then both first and second derivatives of V vanish at </> , ^'(^o) = V"(0o) = 0, and 
the potential becomes extremely flat in the radial direction, see Fig. Jlj]. We note 
that individually none of the terms in eq. ([|) could have driven a successful inflation 
at VEVs lower than Mp = 2.4 x 10 18 GeV. However the combined effect of all the 
terms leads to a successful inflation without the graceful exit problem . 
Around (f> the field is stuck in a plateau with potential energy 



x In what follows we study inflation when eq. (Q) is satisfied. However, this condition is not 
strictly required. Successful inflation is obtained for values of A/Am^ which are slightly larger or 
smaller than 1, (\A 2 — IQmV^I 2 / '4m^) <; f0~ 8 , for detailed discussion see |Ts| ]. 
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The first and second derivatives of the potential vanish, while the third deriva- 
tive does not. Around <fi = 0o one can expand the potential as V{4>) — V(4>o) + 
(l/3!)V"'(0o)(0-0o) 3 , where 

V'"{cf> ) = -j=h 3 rnt. (9) 

Hence, in the range [<p -A<p, (j) +A(j)], where ~ H? ni /V'"((j)o) ~ ((j>o/ M p) > H im, 
the potential is flat along the real direction of the inflaton. Inflation occurs along 
this flat direction. 

If the initial conditions are such that the flat direction starts in the vicinity of (fro 
with <fi « 0, then a sufficiently large number of e-foldings of inflation can be generated. 
Around the saddle point, due to the random fluctuations of the massless field, the 
quantum diffusion is stronger than the classical force, Hi n {/2ii > <p/H in f [jnj, for 



At later times, the evolution is determined by the usual slow roll. The equation of 
motion for the <fi field in the slow-roll approximation is 3H<fi = — (1/2)^ W (0 O )(0 — 



*o) 2 . 



A rough estimate of the number of e-foldings is then given by 

where we have assumed V'(4>) ~ (0 — <fi ) 2 V"' (4> ) (this is justified since ^'(^o) an d 
V'(0o) are both small). We note that the initial displacement from O cannot be 
much smaller than H in {, due to the uncertainty from quantum fluctuations. 

Inflation ends when the slow roll parameters become ~ 1. It turns out that 
1 77 1 ~ 1 gives the dominant condition 



0o 2AhlMl ' 

The total number of e-foldings can be computed as [|TJ: 

L2 \ 1/2 / o X 1/2 



(12) 



x . ( ^ X -( 3m * Y 



(13) 



evaluated after the end of diffusion, see eq. (jTD|), when the slow-roll regime is achieved. 
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Let us now consider adiabatic density perturbations. As in Ref. JE], one finds 

1 Hf nf h\M v 

8b - i ~ -ij A/ CO BE • (14) 

5tt 3m 

In the above expression we have used the slow roll approximation ~ — U'"(0o) (0o — 
(j>) 2 /3H ini , and eq. (O). The number of e-foldings, A/cobej required for the observa- 



tionally relevant perturbations, is > 60 p0| . The exact number depends on the scale 



of inflation and on when the Universe becomes radiation dominated (we note that 
full thermalization is not necessary as it is the relativistic equation of state which 
matters). In our case A/cobe < 60 as we shall see below. 

The spectral tilt of the power spectrum and its running are 

n s = 1 + 277-66-1--^, (15) 

WcOBE 

>hK - ' , (16) 



dink Ac OBE ' 

cf. 0. (We note that e < 1 while r) = -2/A" COB e-) 

It is a remarkable feature of the model that for the weak-scale supersymmetry 
and for the correct value of the Yukawa coupling, namely, 

~ 100 GeV - 10 TeV , h 3 ~ 10" 12 , (17) 

the flat direction N 3 H M L 3 leads to a successful low scale inflation near (f> ~ (10 14 — 10 15 ) GeV 
Mp, with 

V ~ 10 32 - 10 36 GeV 4 , tf inf ~ 10 MeV - lGeV , 

K ~ 10 3 , T max ~ 10 8 - 10 9 GeV . (18) 

The total number of e-foldings driven by the slow roll inflation, A/" e ~ 10 3 , is more 
than sufficient to produce a patch of the Universe with no dangerous relics. Those 
domains that are initially closer to 0o enter self-reproduction in eternal inflation. 
Since the inflaton, N3H U L3, couples directly to MSSM particles, after inflation the 
field oscillates and decays to relativistic MSSM degrees of freedom. The highest 
temperature during reheating is T max ~ V 1 ^. This temperature determines the 
total number of e-foldings required for the relevant perturbations to leave the Hubble 
radius during inflation; in our case it is roughly Acobe ~ 50. 

An interesting observation is that the VEV is smaller than the Planck scale, 
therefore, the potential is free of any trans- Planckian corrections, and the supergravity 
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effects are negligible since 3> H m f . Hence, one can make more certain predictions 
than in the case of the high-scale inflation. 

Despite the low scale of inflation, the flat direction can generate density pertur- 
bations of the correct size for the parameters listed above. Indeed, from eqs. flnyT5| ), 
and ([H])> we obtain: 

dft 

c^~l(r 5 , n s ~0.92, — V 0.002. (19) 

am A; 

The spectral tilt and the running agree with the current WMAP 3-years' data within 
la ||. The tensor modes are negligible because of the low scale of inflation. 

We emphasize that the VEV of the flat direction is related to the Yukawa coupling 
that can generate the Dirac neutrino mass ~ 0.1 eV. The scale of the neutrino mass 
appears to be just right to get the correct amplitude in the CMB perturbations. 

The inflaton has gauge couplings to the electroweak and £/(1)b-l gauge/gaugino 
fields. It therefore induces a VEV-dependent mass ~ g{(f>) for these fields (g denotes 
a typical gauge coupling). After the end of inflation, starts oscillating around 
the global minimum at the origin with a frequency ~ 10 3 iJi n f, see eq. (|18|) . 
When the inflaton passes through the minimum, {(f)) = 0, the induced mass under- 
goes non-adiabatic time variation. This results in non-perturbative particle produc- 



tion ||21|| . As the inflaton VEV is rolling back to its maximum value 0o, the mass 
of the gauge/gaugino quanta increases again. Because of their large couplings they 
quickly decay to the fields which are not coupled to the inflaton, hence massless, 
notably the down- type (s) quarks. This is a very efficient process as a result of which 
the inflaton decays to relativistic particles within few Hubble times after the end of 
inflation (for more details see [0). A thermal bath of MSSM particles is eventually 
formed with a temperature T r h ~ 10 6 GeV (for details of thermalization in SUSY, 
see [22, |23|, 0]). This is sufficiently large to produce cold dark matter in the form 
of thermal neutralinos [E5] . The temperature is also high enough for the electroweak 



baryogenesis |26fl. On the other hand, the reheat temperature is low enough for the 



dangerous relics, such as gravitinos, not to be produced |22, 23, 27 



The inclusion of a gauged £7(1)b-l is important for successful reheating. In the 
absence of it the RH (s)neutrinos would be gauge singlet. Then the energy density 
in the A3 component of the inflaton could only decay via the tiny Yukawa coupling 
fi3, happening at a rate Tn ~ (/i|/87r)m^ ~ 10~ 22 GeV. The oscillations of A3 would 
dominate the universe, and their very late decay (though before BBN, as can be seen 
from r^v) would dilute the generated baryon asymmetry and dark matter. However, 
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the presence of a {7(1)b-l ensures that the energy density in the inflaton completely 
decays through interactions, thus efficiently 

Some comments are in order. In our discussion so far we have included only 
one neutrino and neglected the dynamical role of other two other generations of 



right handed neutrinos, which result in multiple flat directions ||28|| . Their Yukawa 
couplings are even smaller, h 12 <C 10~ 12 . We have checked that the presence of these 
flat directions does not affect the dynamics and the temperature anisotropy. 

We have neglected so far the running of the soft supersymmetry breaking parame- 
ters m.0 and A. In gravity-mediated supersymmetry breaking models these terms run 
logarithmically (as well as the Yukawa coupling h 3 ). Let us focus on the contribution 
of U (1)b-l part. The B — L charges of N, L and H u are +1, —1,0 respectively. Due 
to the smallness of h$ ~ 10~ 12 , the B — L gauge interactions dominate the running. 
Their contribution to the one-loop renormalization group equations of rr&, A, h 3 



closely follows those from the U(l)y piece of the SM [29 



dml 



dA 
^ dfi 
dh?. 



12vr 2 



dfi 



2 I ~ 1 2 i 



167T 



9b—l rh + 



#B_L h 3 + 



(20) 



where \i is the scale, and ^b-l, m are the gauge coupling and gaugino mass of 
the {/(I)b-l respectively. The denotes the contributions from the SM gauge 
interactions, which are similar in size. The running implies VEV dependence of 
m 2 , A, h 3 . This should be taken into account when we determine the flatness 
condition from eq. (^). 

Starting from the grand unified theory (GUT) scale Mqut ~ 2 x 10 16 GeV, we 
have 



A{4>) 



m r 



Ac 



h° 3 



1 + K x In 



1 + K 2 In 



1 + K 3 In 



M 2 



/If 2 
GUT 



/If 2 



(21) 
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Here wig, Aq, h® are the boundary values of parameters at Mqut and 

12vr 2 V m 

*--^Bf) 

^» « -1^3 + - » ( 22 ) 

with mu2 being the gaugino masses at the GUT scale. The again denotes 
contributions of similar size from the SM gauge interactions. 

The condition for the existence of a point 0o such that V'{4>o) = ^"{4>o) = 0> 



and its VEV follow from the analysis of [18 

A(4> ) = 4m^(0 o ) (l + Ki- ^K 2 + l^a) V2 



(l + Ik x - \k,) . (23) 



h° 3 V ' 2 1 4 

Note that these expressions are the same as those given in (|7]||) except for corrections 
due to running. It can be seen from Eq. fl22|) that Ki, K 2 , K 3 ~ 0(1O~ 2 ). Hence 
the corrections are indeed very small. 

However, the important point is that the condition for the flatness of the potential 
involves Tn,$(<f>o), A(<f) ) rather than the boundary values m , A Q . One can then use 
renormalization group equations to translate this into a (slightly different) relation 
between mo and Aq. Therefore quantum corrections do not spoil the flatness of the 
potential; they just slightly modify the flatness condition. In consequence, the rest 
of our analyses remains unchanged: the number of e-foldings ([13]), the amplitude of 
perturbations (0), etc. One should just note that the values of parameters at 0o are 
used in respective expressions 2 . 

Finally the issue of initial condition could be addressed by triggering early bouts 
of inflation as has been often argued, see e.g. Refs. []30], with a last phase 
driven by the A-term inflaton. Further note that near the saddle point of an A-term 
inflation, there is a self-reproduction eternal inflation regime, as we discussed above 
(see eq. (0) and the discussions after eq. (|i~8"D ). The eternal inflation regime along 
with a prior phase of a false vacuum inflation ameliorates the initial condition problem 
considerably for a cosmologically flat direction as we have discussed in Ref. 



2 Note that inflation takes place in an interval A</> ~ lO~ 8 0o around <f>o, see Eq. (|l2|). Therefore 
logarithmic running of parameters in this interval is totally negligible. The only relevant running 
is that from the GUT scale down to 6n ~ 10 14 GeV. 
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To summarize, in the MSSM with right-handed neutrinos the flat direction pa- 
rameterized by the VEV of the right-handed sneutrino, the slepton, and the Higgs 
field can serve as the inflaton. To produce the correct density perturbations, the 
neutrino Yukawa coupling should be small, and the scale of the corresponding Dirac 
masses turns out to be in general agreement with the known neutrino masses. The 
salient feature of the present model is that the inflaton is not an additional ad hoc 
ingredient; it is firmly rooted in particle physics, and its properties can be inferred 
from the upcoming experiments, including the LHC. 
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